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In the context of f{R) modified gravity theories, we study the Kerr-Newman black-hole solutions. 
We study non-zero constant scalar curvature solutions and discuss the metric tensor that satisfies 
the modified field equations. We determine that, in absence of a cosmological constant, the black 
holes existence is determined by the sign of a parameter dependent of the mass, the charge, the spin 
and the scalar curvature. We obtain that for negative values of the curvature, the extremal black 
hole is no longer given by a spin parameter Omax = M (as is the case in General Relativity), but by 
amax < M, and that for positive values of the curvature there are two kinds of extremal black holes: 
the usual one, that occurs for amax > M, and the extreme marginal one, where the exterior (but 
not interior) black hole's horizon vanishes provided that a < a^nin- Thermodynamics for this kind 
of black holes is then studied, as well as their local and global stability. Finally we study different 
f[R) models and see how these properties manifest for their parameters phase space. 



PACS numbers: 98.80.-k, 04.50.-|-h 
I. INTRODUCTION 



General Relativity (GR) has been the most successful 
gravitational theory of the last century, fully accepted 
as a theory that describes the macroscopic geometrical 
properties of spacetime. For an isotropic and homoge- 
neous geometry, GR leads to Friedmann equations which 
describe in an appropriate way the cosmological evolu- 
tion with radiation and then matter dominated epochs. 
Nevertheless, the development of observational cosmol- 
ogy in the last decades with experiments of increasing 
precision like supernovae observations [T], has revealed 
that the Universe is in a stage of accelerated expansion. 
GR provided with usual matter sources is not able to 
explain this phenomenon. Moreover, GR does not ac- 
count either for the cosmological era known as inflation 
[2], believed to have taken place before the radiation 
stage and that could alleviate some problems of stan- 
dard cosmology like the horizon and the flatness prob- 
lem [3 . In addition, GR with usual baryonic matter 
cannot explain the observed matter density determined 
by fitting the standard ACDM model to the WMAP7 
data (Wilkinson Microwave Anisotropy Probe results for 
7 years of observations) [3] , the latest measurements from 
the BAG (Baryon Acoustic Oscillations) in the distribu- 
tion of galaxies and the Hubble constant (i?o) mea- 
surement [5]. Thus, GR requires the introduction of an 
extra component called dark matter (DM), that accounts 
for about 20% of the energy content of our Universe. Al- 
though there are many possible origins for this compo- 
nent [7], DM is usually assumed to be in the form of 
thermal relics that naturally freeze-out with the right 
abundance in many extensions of the standard model of 
particles [8] . Future experiments will be able to discrim- 
inate among the large number of candidates and model. 



such as direct and indirect detection designed explicitly 
for their search [S], or even at high energy colliders, where 
they could be produced [T^ . 

A more puzzling problem is associated to the present 
accelerated expansion of the Universe. There are also 
a large amount of different explanations. One of them, 
assuming the validity of GR, postulates the existence of 
an extra cosmic fluid, the dark energy (DE), whose state 
equation p = wdeP (where p and p are the pressure and 
the energy density of the fluid) demand wde < ^1/3 in 
order to provide an accelerated cosmic expansion (llj . 
The cosmological constant is the simplest model of DE, 
corresponding to an equation of state = — 1. How- 
ever, if we assume that the cosmological constant repre- 
sents the quantum vacuum energy, its value seems to be 
many orders of magnitude bigger than the observed one 

In addition, there is also the problem that, as this cos- 
mological constant cannot account for an inflationary pe- 
riod, a slow rolling scalar field, the inflaton, has to be in- 
troduced by hand. Nevertheless, other explanations for 
the mentioned acceleration may be provided by theories 
that modify GR by considering actions different from the 
Einstein-Hilbert one [T3]. Examples are Lovelock the- 
ories, free of ghosts and whose field equations contain 
second derivatives of the metric at most; string theory 
inspired models, that include a Gauss-Bonnet term in 
the Lagrangian; scalar-tensor theories like Brans-Dicke 
one, in which gravitational interaction is mediated by 
both a scalar field and GR tensor field; or the so called 
f{R) theories, in which our work will be focused. In this 
work we shall restrict ourselves to f{R) theories in the 
metric formalism (where the connection depends on the 
metric, so the present fields in the gravitational sector of 
the action come only from the metric tensor) in the Jor- 
dan frame. In this frame, the gravitational Lagrangian is 
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given by R+f{R), where f{R) is an arbitrary function of 
the scalar curvature R, and Einstein's equations usually 
become of fourth order on the metric derivatives. 

f{R) theories were proved (see [Hj among others) to 
be able to mimic the whole cosmological history, from 
inflation to the actual accelerated expansion era. Di- 
verse applications of these theories on gravitation and 
cosmology have been also widely studied [K], as well as 
multiple ways to observationally and experimentally dis- 
tinguish them from GR. Concerning local tests of gravity 
and other cosmological constraints, see [T5] . 

The study of alternative gravitational theories to GR 
requires to confirm or discard their validity by obtaining 
solutions that can describe correctly, e.g., the cosmolog- 
ical evolution, the growth factor of cosmological pertur- 
bations and the existence of GR-predicted astrophysical 
objects such as black holes (BH). It is a well-known fact 
that, by choosing an appropriate function, f{R) theories 
can mimic any cosmological evolution and, in particular, 
the one described by the ACDM model [T7]. In fact, some 
modified gravity theories present the so called degeneracy 
problem: from large scale observations (la type super- 
nova, BAO, or the cosmic microwave background) which 
depend uniquely on the evolution history of the Universe, 
the nature and the origin of DE cannot be determined 
due to the fact that identical evolutions can be explained 
by a diverse number of theories. However, it has been 
proved [18] that when scalar cosmological perturbations 
are studied, f{R) theories, even mimicking the standard 
cosmological expansion, provide a different matter power 
spectrum from that predicted by the ACDM model [TO] . 
Therefore, it is interesting to study the properties of BH 
in this kind of theories, since some of their known features 
might be either exclusive of Einstein's gravity or intrin- 
sic features of any covariant gravitational theory. On the 
other hand, obtained results could provide a method to 
discard models that disagree with expected physical re- 
sults. In this sense research of BH thermodynamics may 
shed some light about the viability of alternative gravity 
theories since local and global stability regions, and con- 
sequently the existence itself of BH, depend on the values 
of the parameters of the model under consideration. 

BH properties have been widely studied in other mod- 
ified gravity theories: for instance [20l |2T] studied BH 
in Einstein's theory with a Gauss-Bonnet term and a 
cosmological constant. Gauss-Bonnet and/or quadratic 
Riemann interaction terms are studied in |22) . where is 
found that for a negative curvature of the horizon, phase 
transitions might occur. BH in Lovelock gravitational 
theories were studied in |23) . where the corresponding 
entropy was calculated. Other recent works have studied 
[24] BH in the context of Hoi^ava-Lifshitz gravity as well. 
Previous works concerning BH in f{R) theories proved 
that for a Lagrangian R + aR^ the only spherically sym- 
metric solution is Schwarzschild's one provided that one 
works in the Einstein's frame. Again in Einstein's frame, 
[25j proposed uniqueness theorems for spherically sym- 
metric solutions with an arbitrary number of dimensions 



(see for additional results). Spherical solution with 
sources were also studied in [27] whereas |28| developed a 
new covariant formalism to treat spherically symmetric 
spacetimes claiming that Schwarzschild solution is not 
a unique static spherically symmetric solution. Spheri- 
cally symmetric /(i?)-Maxwell and /(i?)-Yang-Mills BH 
were studied on [29, , confirming the existence of numeri- 
cal asymptotic solution for the second ones. Concerning 
axially symmetric solutions, authors in [30] showed that 
these solutions can be derived by generalizing Newman 
and Janis method to f{R) theories. An scalar-tensor 
approach is used in [5T] to show that Kerr BH are unsta- 
ble in a subset of f{R) models because of the superra- 
diant instability. In 32 the entropy of BH is calculated 
in the Palatini formalism by using the Noether charge 
approach. Anti de Sitter (AdS) BH have been studied 
[33] in f{R) models using the Euclidean action method 
(see, e.g., [MIES]) to determine different thermodynamic 
quantities. In [35], the entropy of Schwarzschild-de Sit- 
ter [SdS) BH is calculated in vacuum for certain cos- 
mologically viable models, and their stability discussed. 
In |37| it was proved, in an arbitrary number of dimen- 
sions, that the only static spherically symmetric solu- 
tion -up to second order in perturbations- for a massive 
BH in f{R) theories was that of Schwarzschild-(A)dS'. 
In that same investigation, a thermodynamic analysis of 
Schwarzschild-(A)dS' BH was performed for various f{R) 
models, and it was shown the relation between cosmolog- 
ical and thermodynamic viability. 

This work is organized in the following way: first, some 
general results of f{R) theories in the metric formalism 
are shown in Section II together with the widely accepted 
cosmological viability conditions of f{R) theories. The 
third section is devoted to the study of the axisymmet- 
ric, stationary vacuum solution that describes a massive 
BH with electric charge and angular momentum in these 
theories. In Section IV we study the thermodynamical 
properties of the obtained solutions, whilst fifth section 
analyzes graphically the results of the two previous sec- 
tions for certain f{R) models. Finally, we present the 
conclusions obtained from this work in Section VI. 



II. GENERAL RESULTS 

In order to study the possible solutions obtained from 
any f{R) theory, we start from the action: 

S — Sg -\- Sm J (1) 
where Sg is the gravitational action: 

Sg = ^Jd'x^\{R + f{R)), (2) 

with G = Newton's constant (where Alp is Planck's 

mass), g is the determinant of the metric g^jy (/i, = 
0, 1, 2, 3), R is the scalar curvature of the spacetime and 
f{R) is the function that defines the considered theory. 
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From the matter action term Sm, we define the energy 
momentum tensor as: 

2 6Sm 



(3) 



By performing variations of ([T]) with respect to the met- 
ric tensor, we obtain that the field equations in metric 
formahsm are: 

i?^,(l + /'(i?))-^.g^, [R + f{R)) 

+ (V^V,-5^,n)/'(i?)+87rGT^, = 0, (4) 

with R^i, the Ricci tensor, □ = V/3 (where V is the 
covariant derivative) and f'{R) = df{R)/dR. Taking the 
trace of this equation yields: 

i?(l + /'(i?))-2(i? + /(i?)) -3n/'(i?) + 87rGT =0, 

(5) 

where T — T''^. It is interesting to stress that, unlike in 
GR, vacuum solutions (T — 0) do not necessarily imply 
a null curvature R = 0. From equation Q we obtain the 
condition for vacuum constant scalar curvature R ~ Rq 
solutions: 

(1 + /'(i?o)) - 1 9^. (Ro + f{Ro)) = . (6) 
and the Ricci tensor becomes proportional to the metric: 
^0 + .f{Ro^ 



2(l + /'(i?o)) 



(7) 



with l + f'{Ro) 7^ 0. On the other hand, taking the trace 
on previous equation we obtain: 



Ro (1 + /'(i?o)) - 2 (i?o + /(i?o)) = , 
and therefore 

2/(i?o) 



Ro 



f'iRo) - 1 



(8) 



(9) 



Instead, if f"{R) > 0, a counter reaction mech- 
anism operates to compensate this R growth and 
stabilize the system. 

2.1 + f'{R) > 0. This conditions ensures that the ef- 
fective gravitational constant is positive, as it can 
be checked from the previous definition of Gef / . It 
can also be seen from a quantum point of view as 
the condition that avoids the graviton from becom- 
ing a ghost ^40) . 

3. f'{R) < 0. Keeping in mind the strong restrictions 
of Big Bang nucleosynthesis and cosmic microwave 
background, this condition ensures GR behavior to 
be recovered at early times, that is, f{R)/R — > 
and f'{R) — > as i? — > oo. Conditions 1 and 2 
together demand f{R) to be a monotone increasing 
function between the values — 1 < f'{R) < 0. 

4. f'{R) must be small in recent epochs. This con- 
dition is mandatory in order to satisfy imposed re- 
strictions by local (solar and galactic) gravity tests. 
As the analysis done in [Jl] indicates, the value 
of |/'(i?)| must not be bigger than 10~^ (although 
there is still some controversy about this). This 
is not a needed requirement if the only goal is to 
obtain a model that explains cosmic acceleration. 



III. KERR-NEWMAN BLACK HOLES IN f{R) 
THEORIES 



Since we are looking for constant curvature Rq vacuum 
solutions for fields generated by massive charged objects, 
the appropriate action (in G = c = h = kB — ^ units) is: 



167r 



fiR)-F^,F^n, (10) 



A. Viability conditions of f{R) theories 

The basic conditions and restrictions [38j that are usu- 
ally imposed to f{R) theories to provide consistent both 
gravitational and cosmological evolutions are: 

1- f"{R) > for i? > f'{R). This is the stability re- 
quirement for a high curvature classical regime |39j 
and that of the existence of a matter dominated era 
in cosmological evolution. A simple physical inter- 
pretation can be given to this condition: if an effec- 
tive gravitational constant Ge// = G/(l + /'(i?)) is 
defined, then the sign of its variation with respect 
to R, dGef f I di?, is uniquely determined by the sign 
of f"{R), so in case f"{R) < 0, Ge// would grow as 
R does, because R generates more and more curva- 
ture itself. This mechanism would destabilize the 
theory, as it wouldn't have a fundamental state be- 
cause any small curvature would grow to infinite. 



where Fj^^^ — d^^ Ai, ~ d^A^ and the electromagnetic 
potential. This action leads to the field equations: 

i?^,. (1 + /'(i?o)) - I .9m. {Ro + f{Ro)) 



-2 l^F^^Fl - -g^^F^pF'^P j = . (11) 

At this stage, it is worth stressing that if wc take the 
trace of the previous equation, ([8| is recovered due to 
the fact that i^'f, = 0. 

The axisymmetric, stationary and constant curvature 
Rq solution that describes a BH with mass, electric 
charge and angular momentum was found by Garter 
and published for the first time in 1973 [12]. In Boyer- 
Lindquist coordinates, the metric describing with no co- 
ordinate singularities the spacetime exterior to the BH 
and interior to the cosmological horizon (provided it ex- 
ists, as will be studied below), takes the form: 
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Ae sin 



dt 



Ar f dt 



d(j)\ 



, „ — a sm" w — 



(12) 



r 



with: 

: 



a-) 



Rq 

I2 



2Mr 



(l + /'(i?o))' 



+ cos^ i 



12 



1 



T2 



(13) 



where M, a and g denote the mass, spin and electric 
charge parameters respectively. Notice that, unhke in 
the GR case, the contribution of the charge of the BH to 

the metric is corrected by a (1 + /'(i?o)) factor. This 
feature was aheady obtained for Reissner-Nordstrom BH 
in [31]. 

On the other hand, the required potential vector and 
electromagnetic field tensor in equation ( 11 ) solutions for 
metric ( 12 ) are respectively: 



F = - 



qr f dt . 2qM 
— — asm 6/ — 

q (r^ — cos^ ( dt 

7 Is 

2qr a cos 6 sin 9 , „ 
— de A 



— a sm 



A dr 



a-T^ — [r + a ) 



(14) 

To lighten notation, from now on we will use = 
q^ I (l + /'(i?o)) to refer to the electric charge param- 
eter of the BH. 



The nature of coordinates in (12 1 can be seen by con- 
sidering the M 0, Q ^ 0, i?o 
metric. Thus ( 12 ) becomes: 



limits on this 



ds^ = -dt'- 



P 



dr^ +p2d6|2 + {r^ + a^)sm^ 



(15) 



i.e., Minkowski spacetime in spacial coordinates 



(r, 6*, 0). It is not obvious from (15) that what one has 



is Minkowski spacetime, this is because Boyer-Lindquist 
coordinates need to be "untwisted" via cartesian coordi- 
nates to confirm that what we have is actually an empty 
spacetime: 



sm U cos ( 



r cos 9 



sm tf sm ( 



(16) 



with r>0, 0<6'<7ryO<(/>< 27r. Nevertheless, one 
must keep in mind that, when M 7^ 0, Q 7^ and i?o 7^ 0, 
the simplest interpretation given to these coordinates is 
not completely appropriate due to the distortion of the 
empty spacetime that the presence of the BH induces. 
On the other hand, if we do M -> 0, Q — > 0, a — !> 0, we 
obtain a constant curvature spacetime metric: 



ds^ 



1 



Rq 
12 



dt' 



1 



1 - 



Rpr' 
12 



dr' 



(17) 



that corresponds to either dS or AdS spacetime de- 
pending on the sign of Rq. It is also easy to verify that 
when a — 0, Q — >■ 0, Schwarzschild-(A)(iiS' BH is recov- 
ered. 



A. Singularities 

We will study now the singularities of these BH. Cal- 
culating R^'"^'' R^uap, only p — Q happens to be an in- 
trinsic singularity, and considering the definition of p in 
(13), such singularity is given by: 



and 



7r/2. 



(18) 



Keeping in mind that we are working with Boyer- 
Lindquist coordinates, the set of points given by r = 
and 9 — TT /2 represent a ring in the equatorial plane of 
radius a centered on the rotation axis of the BH, just as 
it happens in Kerr BH [43j. 



B. Horizons 

It is also interesting to study the horizon structure of 
these BH: according to the horizon definition g^'^ ~ 0, 
they are found as the roots of the equation — 0, that 
is: 



a'- 



12 

Rp, 



24 M 

Rn 



■'-R-J^ 



Q') = 0, 



(19) 



fourth order equation that can be rewritten as: 

{r - r^){r - rint){r - rext)ir - rcosm) =0, (20) 

where r^ is always a negative solution with no physi- 
cal meaning, rint and rext are the interior and exterior 
horizon respectively, and rcosm represents - provided it 
arises, as will be seen later - the cosmological event hori- 
zon for observers between r^xt 
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Figure 1. Graphics sliowing positions of liorizons as solutions of tiic equation A,- = 0. On tiie left panel (Rg < 0) we show the cases h > 
(I, BH with well-defined horizons, dashed with dots), h = (II, extremal BH, continuous line) and h < (III, naked singularity, dashed). 
On the right panel {Rq > 0) we represent the cases /i < (I, BH with well-defined horizons, dashed with dots), h = (II, extremal BH 
and III, extremal marginal BH, continuous line), and h > (IV, naked singularity and V, naked marginal singularity, dashed). 



divides the region that the observer could see from the 
region he could never see if he waited long enough time. 



Using L. Ferrari's method |44j to solve quartic equations, 
the existence of real solutions for this equation is given 
by a factor h let us name it horizon parameter. 



h = 



4 A fio 2 
— 1 a 

Ro V 12 



H 3 



-4 a' 




- 18 



(21) 



For a negative scalar curvature Rq, three options may 
be considered: i) h > 0: there are only two real solutions, 
rint and rea;t, lacking this configuration a cosmological 
horizon, as it is expected for an AdS like Universe, ii) 
h = 0: there is only a degenerated root, particular case 
of an extremal BH, whose interior and exterior horizons 
have merged into one single horizon with a null surface 
gravity k (that will be defined in section IV). in) h < 0: 



it is found that there is no real solution to (21), which 
translates into an absence of horizons that leads to a 
naked singularity. 

For a positive curvature Rq , there are also several con- 
figurations depending on the value of h: i) h < 0: both 
^mtj Text and rcosm are positive and real, thus the BH 
possesses a well-defined horizon structure in an Universe 
with a cosmological horizon, ii) h ~ Q: two different 
cases may be described, either rmt and rext become de- 
generated solutions, or rg^jj and Tcosm do so. The first 
case represents an extremal BH, described before. The 
second one can be understood as the cosmological limit 
for which a BH preserves its exterior horizon without be- 
ing "torn apart" due to the relative recession speed be- 
tween two radially separated points induced by the cos- 



mic expansion in an Universe described by a constant 
positive curvature; this case is known as marginal naked 
singularity. Hi) h > 0: There is only one positive root, 
that may be either rint or rcosm- In the first case, the 
mass of the BH has exceeded the limit imposed by the 
cosmology (just described for h — 0), and there are nei- 
ther exterior nor cosmological horizon. This situation 
just leaves the interior horizon to cover the singularity 
{marginal naked singularity case. On the contrary, if the 
root corresponds to r^osnn this time there is a naked sin- 
gularity with a cosmological horizon. In Figure[l]we have 
shown the zeros of the fourth order polynomial A^: it is 
remarkable that, from a certain positive value of the cur- 
vature Rq"* onward, the h factor goes to zero for two 
values of a, i.e., apart from the usual Omax for which 
the BH turns extremal, there is now a spin lower bound 
imim below which the BH turns into a marginal extremal 
BH, as we discussed before. Therefore 



h{a„,ax, M, \Ro\ >0, g) =0 

{M,\Ro\>0,Q), (22) 
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-0.5 0.0 

Figure 2. The shaded regions, delimited by the upper amax and lower amin curves, represent the values of a/AI for which the existence 
of BH is possible once Rg M"^ value is fixed. It is shown for Q / M = and Q / M = 0.75 on the left and right panels repectively. Note 
that Rq has dimensions of [length]"^ when normalizing. 



(M, i?o > i?o'"* > 0, Q) . (23) 



Due to the excessive length of the equations that describe 
the behavior of amax and Omm, we prefer not to display 
them here. Instead, in Figure [2] we show -for certain 
values of the electric charge Q parameter Q - the range of 
values of the spin a parameter for which BH are allowed 
taking into account. To do so, the corresponding Rq 
value is determined by the parameters defining each f{R) 



model, as can be seen from equation (|8|) 



C. Stationary Limit Surfaces 

Another interesting feature of Kerr-Newman BH are 
Stationary Limit Surfaces (SLS), given by gu — 0. For 
Boyer-Lindquist coordinates, this condition translates 
into: 



p2 ■E? 



0, 



(24) 



that leads to the fourth order equation 



12 



24 M 
Rq 



2 2 , 

a cos I 



12 



SVC? 9 



12 

Rq 



(25) 



which can be rewritten as: 

(r - rs-){r - rstnt){r - rsext){r - rscosm) = 0.(26) 

From this equation it follows that each horizon has one 
"associated" SLS. Both hypersurfaces coincide at 9 — 
0,7r as seen when comparing (25 1 with equation (19). A 



complete scheme of BH horizons and SLS structure is 
shown in Figure [3] for both signs of Rq. 



IV. BLACK HOLE THERMODYNAMICS 

From now on, we will focus our study on BH with a 
well-defined horizon structure and only for Rq negative 
values of. This last choice is motivated by the problems 



arising when normalizing the temporal Killing ^ = dt in 
positive curvature spacetimes. This problem is more ex- 
tensively examined in f?S] . i?o < choice will allow us to 
define the thermodynamical quantities of the BH. The ro- 
tation Killing field ■0 = is uniquely determined by the 
condition that their orbits should be closed curves with 
a length parameter equal to 27r. Nevertheless, there is 
not an adequate criterion to normalize the Killing vector 
^ = £?t in the dS {Rq > 0) Universe since multiplicative 
constants can be added to ^ and the obtained Killing 
vectors continue being null on the horizon. In the AdS 
{Rq < 0), the normalization is done without problems by 

1/2 

imposing that the ^ value tends to r (—12 / i?o) as r 
goes to infinity. 

In order to study the different thermodynamical prop- 
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Exterior horizon 
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Interior SLS Interior horizon 




Cosmologioal SLS 



Figure 3. On the left: diagram of a Kerr-Newman BH structure witli negative curvature solution Rg = —0.4 < 0, A/ = 1, a = 0.85 and 
Q = 0.35 (h > 0). Right: BH structure with positive curvature solution Ro = 0.4 > 0, M = 1, a = 0.9 and Q = 0.4 {h < 0). Dotted 
surfaces represent the SLS whereas horizons are shown with continuous lines. The rotation axis of the BH is indicated by the vertical 
arrow. In both types of BH, the region between the exterior SLS ext and its associated exterior horizon r^xt is known as ergoregion. 



erties of Kerr-Newman BH in f{R) theories, we start 
looking for the temperature of the exterior horizon Vext = 
Text {Rq, Q, M). For that pm^pose, we will use the 
Euclidean action method 46 . Performing the change 
t — ^ —IT, a — ;> ia on the metric ( 12 ) we obtain the Eu- 



clidean section, whose non singular metric is positive- 
definite, and time coordinate has now angular character 
around the "axis" r = r^xt ■ Regularity of the metric at 
T = Text requires the identification of points: 



{r, r, 0, 



{T + ip,r,0,^ + ipnH), (27) 



where /3, which represents the period of the imaginary 
time on the Euclidean section, it's the inverse Hawking 
temperature: 



13 



47r (r 



Text 



Ro Ro Tl^t (a2 + Q^) 



12 



1 



(28) 



and Qh is the angular velocity of the rotating horizon, 
which is the same for all the horizon: 



H 



(29) 



Considering that the event horizon is also a Killing hori- 
zon of the Killing vector x = C + ^E'lp (where, as was 
said before, ^ = dt and = are the vectors that 
asymptotically represent time translations and rotations 
respectively), fin could be also obtained demanding x to 
be a null vector on the horizon: 



= 



(30) 



BH horizon temperature could have also been obtained 
through Killing vectors, as is explained in [47j where tem- 
perature is defined as follows: 



K 

4^' 



with K the surface gravity defined by: 



KXiy- 



(31) 



(32) 



It can be verified that k is the same at any horizon point 
and consequently = Te as obtained in 148] . 

Now that we know the expression for the temperature, 
we consider the Euclidean action in order to obtain the 
remaining thermodynamical quantities: 



ASe 



g\{Ro + /(i?o) 



(33) 



with y the integration region. As is described in [53], one 
has to calculate the difference in four- volumes of the two 
metrics, identified by the same imaginary time. Provided 
that the metric is stationary, integration over time simply 
leads to a multiplicative factor f3. On the other hand, 
keeping in mind that Maxwell's equations V ^F^'^ = 
must be satisfied, we can rewrite the third term in the 
integrand as a divergence: 

F^.F^'' ^V,{2F>''' A^), (34) 

and therefore: 

16 TT StT Jgy 



(35) 



8 



where dy = 3^x3"^ is the boundary of the considered re- 
gion, is a 2-sphere whose radius has to be sent to infin- 
ity after the integration, and AV is the volume difference 



between both sohitions (corresponding to the black hole 
metric and that of AdS space identified with the same 
imaginary time). After some calculation, wc obtain: 



P iRo + f{Ro)) 
24 S 













Ro) 



12 a' 

Ro rea 



^ ' ext 



(36) 



r 



where $e is the electric potential of the horizon as seen 
from infinity: 

$e = A.x'-l^^ - A.x'l^^^^^ = -ff-^ , (37) 

and Q is the physical electric charge of the BH, obtained 
integrating the fiux of the electromagnetic field tensor at 
infinity, which happens to be: 



Q = ^ 



(38) 



We shall remember that these calculations involve the 
vector potential and the electromagnetic field tensor 

, . 1/2 

given in ( 14 1, and that's why the factor (1 -I- /'(i?o)) 
does not appear here. Further analysis of the action re- 
veals that it goes singular for /i = 0, as could be expected 
from extremal BH, whose temperature Te = makes the 
/? factor diverge. Since thermodynamical potentials are 
obtained by dividing the action by the (3 factor, they still 
remain well defined at — 0. 



It can also be seen that the action ( 36 1 diverges in the 
2 — 12 /Rq, which implies 5 = 0. This singu- 



limit a 

lar case ^49] is further explored in fSO^ and implies that 
a 3-dimensional static closed Universe at infinity would 
rotate with the speed of light. Thus, in order to avoid 
all these problematic issues, let us assume from now on 
that: 



S := 1 



12 



> 0. 



(39) 



The previous expression will turn out also to be a re- 
quired condition to ensure a positive area and entropy of 
the BH, as will be seen below. 



By using the expression ( 36 ) we can immediately ob- 



tain Helmholtz free energy F, defined by: 



F : 



ASe 



(40) 



where the term fljj J comes from the required Legen- 
dre transformation to fix angular momentum, being J 
the angular momentum of the BH and Qh the angular 
velocity of the horizon computed before in (29). To cal- 



culate J we need first the physical mass associated to the 
BH, which can be calculated from: 



resulting on an angular momentum: 



J = 



aM 



aM , 



r{Ro)) , 



(41) 



(42) 



where we have used equation ([S]) on Ai to make the sub- 
stitution: 2 {Ro + /(i?o)) / i?o = 1 + /'(-Ro)- Using again 
the relation Q"^ — j (1 + f'{Ro)), and equation (19 1 to 
express M as a function of r^xt, we obtain: 



F=(l + /'(i?o)) 



[36 



12 rf 



Ro + a'^ (36 



Ro)] 



24:rext^ 



(43) 



If the condition l + /'(_Ro) > is required to hold in order 
to obtain positive values of the mass, by analyzing the 
numerator of i^, we find F > for values of the horizon 
(with an associated mass M'™'* through 



below r ''™** 



equation (19)), and F < for larger values. 

Using the appropriate thermodynamical relations |52j , 
we can derive the entropy S of the BH, which reads: 



S = l3{M-nHJ)-ASE 

= ii+nRo))^-^^k^. 



(44) 



If we compute now the area Ah of the exterior horizon 



Text, which can be calculated from the metric (12) doing 
"f" = fext and t constants, we obtain: 



dsi, = 



A, 



A 



H 



det g(2) de 



-(rext + a ) — 



4^(^?. 



(45) 



(46) 



Therefore one sees straightforwardly that the entropy 
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(44 1 can be expressed as: 



5=(l + /'(i?o)) 



A 



H 



(47) 



consequently l+f'{Ro) > is also a mandatory condition 
to obtain a positive entropy [ST], as we supposed above. 

Once the temperature T and the entropy S of the BH 
are known, we can take a step further and study the heat 



capacity at constant scalar curvature Rq and at fixed spin 
a and charge Q parameters. From the definition: 



C = T 



dS 
df 



(48) 



Ro,a,Q 



we obtain the expression: 



C = (1 + /'(i?o)) - 



27r {a^ + r^t) W (12 + r^t i?o) + 3 (4 



4r'^ 



Ro)] 



E [-36 Q2 ^2^^ + ^4 (_i2 + i?o) + 3 r^^ 



r 



(4- 



Ro)-ia' (3Q2 + i2r2^^_2rLt i?o)] 



(49) 



Provided that the condition (39) holds, it seems interest 



ing to find out for which values of Rq, a, Q and M the 
denominator of the thermal capacity goes to zero, i.e., 
the thermal capacity goes through an infinite disconti- 
nuity, which corresponds to a BH phase transition. We 
can distinguish between two kind of BH on this subject 
depending on the values of the a, Q and M parameters 
and scalar curvature Rq: i) fast BH, without phase tran- 
sitions and always positive heat capacity C > 0. ii) slow 
BH, presents two phase transitions for two determined 

values of Text- 



dT 



t) = 0, 



(50) 



f9) 



lehavior of the tem- 



with Tg^j < rg)j.f , being r\^^ a local maximum of the tem- 
perature Tmax, and r^^^^ a local minimum Tmin- BH heat 
capacity proves to be positive (C > 0) for r^^t < i"\xt ^^"^ 
Text > r^exu and negative (C < 0) for rl^^ < r^xt < r^ext- 
Once again two masses and AI^^ can be associated to 
the radii rl^^ and r^^^ via equation 

In Figure|4]we have visualized the 
perature T, the free energy F and the heat capacity C of 
a BH for different values of mass M, with fixed a, Q and 
Rq values [? ] . It can also be seen the range of a and 
Q parameters values that provide slow or fast BH for a 
constant value of the scalar curvature Rq. 

Unlike Schwarzschild-AdS' BH case [34], Kerr- 
Newman-AdS' BH are allowed for any value of the tem- 
perature T, hence stability of each BH will be exclusively 
given by the corresponding values of heat capacity C and 
free energy F as functions of the a, Q, M and Rq param- 
eters, that ultimately define the BH. However, for a set of 
fixed values of a, Q and Rq, the mass parameter must be 
bigger than a minimum M™™ (characterized by T = 0) 
to have BH configuration, otherwise radiation is the only 
possible equilibrium up to such a minimum mass. For 
bigger masses, wc shall distinguish between the fast and 
the slow BH. Fast BH, with bigger values of the spin and 
the electric charge than the slow ones, shows a heat ca- 
pacity always positive and a positive free energy up to 
a, M = M '*™** value, and negative onwards. Thus, this 



BH is unstable against tunneling decay into radiation for 
mass parameter values of M < M For M > M 
free energy becomes negative, therefore smaller than that 
of pure radiation, that will tend to collapse to the BH 
configuration in equilibrium with thermal radiation. 

The second situation, i.e., the slow BH shows a more 
complex thermodynamics, being necessary to distinguish 
between four regions delimited by the mass parameter 
values: M""™ < < M" < M'*™*. For M™" < M < 
and for M" < M < M both the heat capacity 
and the free energy are positive, which means that the 
BH is unstable to decay by tunneling into radiation. If 
< M < the heat capacity becomes negative 

but free energy remains positive, being therefore unstable 
to decay into pure thermal radiation or to larger values 
of mass. Finally, for M > Af'™'* the heat capacity is 
positive whereas the free energy is now negative, thus 
tending pure radiation to tunnel to the BH configuration 
in equilibrium with thermal radiation. 

It is mandatory to say that, although not quantita- 
tively, the thermodynamical behavior of these f{R) BH 
is qualitatively similar to that of GR [S2] . 



V. PARTICULAR EXAMPLES 

In this section we will study some particular f{R) mod- 
els. For each model we will firstly study the range of pa- 
rameters that allows the existence of Kerr-Newman BH. 
Secondly, we will focus on the thermodynamical quanti- 
ties that define BH stability depending again upon the 
model range of parameters [M] . 

For the sake of simplicity, let us introduce the dimen- 
sionless variables: 



r 

M 



a 

M 



M 



->Q, RqAI^^Rq, (51) 



where Ad is the mass parameter, a the spin parameter, 
Q the electric charge parameter and Rq is the scalar cur- 
vature obtained as a solution of equation ([s]). The con- 
sidered models are: 
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Figure 4. For Rq = —0.2, we graphically display temperature (up to the left), heat capacity (up to the right), and the free energy (down 
to the left) of a BH as functions of the mass parameter M for the cases: I) a = 0.5 y Q = : "slow" BH that shows a local maximum 
temperature Tmax and a local minimum temperature Tmin at the points where the heat capacity diverges, taking the latter negative 
values between Tmax y Tmin- H) i ~ 0.965 y Q = : limit case where Tmax and Tmin merge, hence resulting on an inflection point in the 
temperature and an always positive heat capacity. Ill) a = 1.5 y Q = : "fast" BH with both temperature and heat capacity monotone 
growing (always positive too). It can be seen that all the configurations acquire a value F < from a certain value of M onward, given by 
rj"™'. The values of Af™'", with T = and C = 0, correspond to an extremal BH. Down to the right, for Ro = —0.2 aswell, we display 
the regions in which BH behave as "slow" or "fast" BH. 



A. Model I: f{R) = a\Rf 

This model has been widely studied because the aR^ 
term with a > can account for the accelerated ex- 
pansion of the Universe. This model can also explain the 
observed temperature anisotropics observed in the CMB, 
and could become a viable alternative to scalar field in- 
flationary models; reheating after inflation would have its 
origin on the production of particles during the oscilla- 
tion phase of the Ricci scalar [SS]. By using expression 
([9]), we obtain the following scalar curvature solutions: 



where solution leads to a positive curvature, and 
to a negative one. The viability condition 1 + f'{R^) > 
restricts the range of parameters that define this f{R) 
model to different regions depending on what solu- 
tion we choose, R^ or Rq-, for R^ we have Region 1 
{a > 0, /3 > 2} and Region 2 {a < 0, l3 < 1}, and for i?g , 
Region 3 {a < 0, /3 > 2} and Region 4 {a > 0, /3 < 1}. 
In Figure [5] we show the range of the spin parameter a 
for which BH are allowed, depending on the parameters 
a y /3 and for certain values of the charge parameter Q. 
We graphically schematize in Figure [9] the possible ther- 
modynamical configurations as functions of a and /3, for 
those regions in which i?o < 0- 



±1 



(/?-2), 



(52) 
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B. Model II: f{R) = ±li?l"exp 



R 



R 



For this model and by using again expression ([9]), the 
scalar curvature, independently of the chosen sign, be- 



comes: 



/3 



a-2 



(53) 



However, the condition 1 + /'(i?o) > limits the theory 
f{R) to different regions depending on what sign we de- 
cide to work with. If Rq > 0, the theory is restricted to 
Region 1 {a > 2, ^ > 0} and Region 2 {a < 2, /? < 0}, 
for which Rq takes positive values. If Rq < 0, we re- 
strict ourselves to Region 3{a>2, /3<0} and Region 4 
{a < 2, /? > 0}. In Figure [6] we show, just as before, the 
spin parameter a values for which BH can exist in this 
model depending on the a and /3 parameters defining the 
model. We graphically schematize in Figure [9] the pos- 
sible thermodynamical configurations as functions of a 
and /3, for those regions where i?o < 0. 



C. Model III: f{R) = R{hog{aR)f - R 



The associated scalar curvature to this model is: 



Rq = -exp(/3) 
a 



(54) 



In this case the condition 1 + f'{Ro) > restricts us to 
work with Region 1 {a e M, /3 > 0}, where Rq takes pos- 
itive values for a > and negatives for a < 0. In Figure 
[T] we graphically represent spin parameter a values for 
which BH present their complete horizon structure, de- 
pending on the values of the parameters that define the 
model. We graphically schematize in Figure [9] the differ- 
ent possible thermodynamical configurations as functions 
of a and /3, for those regions where i?o < 0. 



D. Modelo IV: f(R) = -a 



R 



1 + 



R 



This model has been proposed [SS] as cosmologically 
viable. For our study we will consider the case n = 1, 
thus having a biparametric theory, as we can define 7 = 
/3/a and then obtain: 



fiR) 



kR 

l + lR 



(55) 



Replacing the latter in ^ we obtain two different values 
for the curvature: 



R^ 



1 - K 



± 



-K, (1 — k) 



(56) 



Keeping in mind that we have to satisfy 1 -I- f'{Ro) > 0, 
K happens to be restricted to values k > 1. On the 
other hand, computation of 1 + /'(Rq ) reveals that Rq 
is only a valid solution for values of k and 7 in the Re- 
gion 1: {k > 1, 7 > 0}, and Rq only in Region 2: 
{k > 1, 7 < 0}, being Rq > and Rq < in their 
respective regions. In Figure [8] we show the range of the 
spin parameter a for which BH are allowed, depending on 
the parameters k and 7 for certain values of the charge 
Q. We graphically schematize in Figure [9] the different 
possible thermodynamical configurations as functions of 
K and 7, for those regions in which _Ro < 0. 



VI. CONCLUSIONS 

In this work we have derived the metric tensor that 
describes a massive, charged, spinning object for f{R) 
gravity in metric formalism. We found that it differs 
from that found by Carter |32] by a multiplicative factor 
in the electric charge and a redefinition of vacuum scalar 
curvature. 

Further study of the metric allowed us to describe the 
different astrophysical objects derived from the presence 
of different horizons: BH, extremal BH, marginal ex- 
tremal BH, naked singularities and naked extremal sin- 
gularities. Focusing on BH and their horizon structure, 
we have seen that these can only exist for values of the 
spin lower than a maximum value amax, and that from 
a certain positive value of the curvature onward, only 
above a minimum value amin- We have also studied the 
thermodynamics of ^dS'-like BH (negative curvature so- 
lutions) by employing the Euclidean action method. It 
is observed that some quantities such as the mass, the 
energy or the entropy of these BH differ from those pre- 
dicted in GR by a multiplicative factor 1 -I- /'{Rq). This 
factor has to be positive in order to assure a positive mass 
and entropy for these kind of BH. 

On the other hand, analysis of the behavior of the 
heat capacity of these BH reveals that we can distinguish 
between two kind of BH: fast and slow, showing the 
latter two phase transitions. We have also investigated 
the stability of the different possible configurations that 
arise from the values of the free energy and the heat 
capacity, observing that qualitatively the situation is 
similar to that described by Kerr-Newman-AdS" BH. 
Finally, we considered four /(i?) models and analyzed 
graphically the previous obtained results. Experimental 
checks to test the validity of a particular f{R) model 
might be done not only by studying astrophysical BH 
stability but also if quantum gravity scale is near TeV: 
since LHC would be producing about one microBH per 
second [57j, stability and thermodynamical properties 
of the produced BH might shed some light about the 
underlying theory of gravity. In this sense, let us remind 
that the relation between BH mass and temperature in 
f{R) theories would differ from that predicted by GR. 
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Figure 5. Model I. Region 1: {a > 0, (5 > 2}, Region 2; {a < 0, /3 < 1}, Region 3: {« < 0, /3 > 2} and Region 4: {a > 0, /3 < 1}. BH 
witli a well defined horizon structure will only exist if they have a spin parameter below the upper surfa<;e amax, and above a second 
surfax;e amin (in case it exists, Regions 1 and 2 for this model) for certain values of a and 
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Figure 6. Model II. Region 1: {a > 2, /3 > 0}, Region 2: {a<2, < 0}, Region 3: {q > 2, /3 < 0} and Region 4: {a < 2, /9 > 0}. 
BH witli a well defined horizon structure will only exist if they have a spin parameter below the upper surfa<;e Umax, and above a second 
surfax;e amin (in case it exists, Regions 1 and 2 for this model) for certain values of a and 
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Figure 7. Model III. Region 1: {a S K, /3 > 0}. BH with a well defined horizon structure will only exist if they have a spin parameter 
below the upper surface amax, and above a second surface a-min (that appears in this case for values of a > 0) for certain values of a and 
/3. 



Model rv, RegiDn 1, Q = 




Figure 8. Model IV. Region 1; {k > 1, 7 > 0} and Region 2: {k > 1, 7 < 0}. BH with a well defined horizon structure will only exist 
if they have a spin parameter below the upper surface amax, and above a second surface amin (in case it exists, only in Region 1 for this 
model) for certain values of k and 7. 
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Figure 9. Thcrmodynamical regions witii negative scalar curvature _R() < of the models I, II, III and IV. For the sake of simplicity, we 
study a BH with the following parameter values: M = 1, a = 0.4 and Q = 0.2. We distinguish between three different regions: i) C < 
and F > 0, in black, ii) C> and F > 0, in gray, iii) C >0 and F < 0, in white. 



